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Abstract 

Sharp upper bounds for the first eigenvalue of the Laplacian on a 
surface of a fixed area are known only in genera zero and one. We 
investigate the genus two case and conjecture that the first eigenvalue is 
maximized on a singular surface which is realized as a double branched 
covering over a sphere. The six ramification points are chosen in such 
a way that this surface has a complex structure of the Bolza surface. 
We prove that our conjecture follows from a lower bound on the first 
eigenvalue of a certain mixed Dirichlet-Neumann boundary value problem 
on a half-disk. The latter can be studied numerically, and we present 
conclusive evidence supporting the conjecture. 
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1 Introduction and main results 



1.1 Upper bounds on the first eigenvalue 



Let M be a closed surface of genus 7 and let g be the Riemannian metric on 
M. Denote by A the Laplace-Beltrami operator on M, and by A x the smallest 
positive eigenvalue of the Laplacian. Let the area Area(M) be fixed. How large 
can Ai be on such a surface? 

Sharp bounds for the first eigenvalue are known only for the sphere ( [HJ , see 
also |5E|), the projective plane ([EZ]), the torus ([Be7|, 03] ), and the Klein bottle 
( |JNPj . jEGJj ). The present paper is concerned with the surface of genus 2. 

Let M be orientable and let II : M — > S 2 be a non-constant holomorphic 
map (or, conformal branched covering) of degree d. It was proved in |YY| that 



AiArea(M) < 8nd . 



(1.1.1) 



Any Riemann surface of genus 7 can be represented as a branched cover over 



S> 2 of degree d 
Therefore, 



where [■] denotes the integer part (see jGunj . pHJ). 

~7 + 3" 



AiArea(M) < 8w 



(1.1.2) 



In general, (|1.1.2|) is not sharp, for example for 7 = 1 (|NJ). 

Let M = V be a surface of genus 7 = 2. Then (ll.l.2|) implies 

AiArea(P) < 16vr. 



(1.1.3) 



The aim of this paper is to show, using a mixture of analytic and numerical tools, 
that ()1.1.3|) is sharp. Main results of this paper were announced (without proofs) 
in jJLNPI section 4]. 



1.2 The Bolza surface 

Let II : V — > § 2 be a branched covering of degree d = 2. The Riemann-Hurwitz 
formula (see |GHj) implies that this cover is ramified at 6 points. We choose 
these points to be the intersections of the round sphere S 2 centered at the origin 
with the coordinate axes in R 3 . The surface V can be realized as 

{(,,«,) EC:^F( ;) : = m |. 
I (z+l)(z + i)J 
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This surface has the conformal structure of the Bolza surface. It has an octahe- 
dral group of holomorphic automorphisms and its symmetry group is the largest 
among surfaces of genus two [11 IKWj . Interestingly enough, the Bolza surface 
appears in some other extremal problems, in particular for systoles (see |K5| ). 

To simplify calculations it is convenient to rotate the equatorial plane by 7r/4. 
The equation of V becomes 

V := (z,w) E C 2 : w 2 = F(z) := z ) m \ J- . (1.2.1) 

The projection II is defined by II : (z,w) — > z. The set of ramification points in 
the complex z plane is R := {0, oo, ie 7 ™/ 4 , ±e 3m ^}. The spherical and complex 
models are related by the stereographic projection; the induced metric in the 
complex plane (which we assume coincides with the equatorial plane of § 2 ) is 

Adzdz/(1 + \z\ 2 ) 2 . (1.2.2) 

Let go be the metric on V which is the pullback of the round metric ()1.2.2|) 
on S 2 . One can see that the metric go has conical singularities at the points 
of ramification. It has curvature +1 everywhere except the branching points. 
Because of the presence of singularities we have to specify what we mean by the 
first positive eigenvalue of the Laplacian on (V,go). We set 

II Vn|| 2 

AiCP,#o):= , mf 



uei^CP.so), u^o, (u,i)=o \\u\\ 2 

where the scalar product (•, •) and the norm || • || are taken in the space L 2 (V, go)- 
The Sobolev space Hl(V,go) of functions supported away from the singularities 
is obtained by the closure of C^(V, g ) := {v G C°°(V, g ) \ U supp v n R = 0} 
with respect to the norm || Vv|| 2 + ||w|| 2 . 

1.3 Main results 

We start with the following 

Conjecture 1.3.1. The equality in (II .1.3(1 is attained for the metric g on V , 
i.e. 

X 1 (V,go)Aie a (P,go) = 16tt. 
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Since (V,go) is a double cover of the standard § 2 , we have 
Area(P, g Q ) = 2Area(§ 2 ) = 8tt. 
Therefore, in order to prove Conjecture 11.3.11 it suffices to show that 

Ai(P,<? ) =Ax(§ 2 ) =2. (1.3.2) 

Unfortunately, we are unable to prove (jl.3.2j) . and therefore establish Con- 
jecture H3III We can however reduce the conjecture to the following spectral 
problem on a quarter-sphere Q C § 2 that can be treated using numerical meth- 
ods. Namely, let, in usual spherical coordinates (4>,9), 

Q = {(<p,9) : < <p < 7r/2, < 9 < tt} . 

We split the boundary dQ into two parts: dQ = d\Q U d 2 Q, where 

d x Q = {(0, 9):\9- ir/2\ < tt/4} U {(tt/2, 9): 0<9< tt/2} , 

d 2 Q = {(0, 9): \9- tt/2| > vr/4} U {(tt/2, 9) : tt/2 < < tt} , 

and consider the spectral boundary value problem for the Laplace-Beltrami op- 
erator on Q: 

-Au = Au on Q , u\ dlQ = , (du/dn)\ d2 Q = . (1.3.3) 

Let Ai denote the first eigenvalue of the problem Q1.3.3j) (which we under- 
stand as usual in the variational sense). 

Conjecture 1.3.4. 

Ai > 2. 

Our main result is 
Theorem 1.3.5. Coniecture \1.3.4\ implies Coniecture \l.3.1l 
Theorem 11.3.51 is proved in sections 2 and 3. 

Extensive numerical calculations (see Section |4j) show that Ai > 2.27 which 
implies Coniecture ll.3.11 The best lower bound we are able to prove is just Ai > 
0.75, which follows from Dirichlet-Neumann bracketing (replace the Dirichlet 
condition by the Neumann one on the arc (0,0), 7r/4 < < 37r/4), see |Ke]). 
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Finally, we note that the spectral problem (jl.3.3|) easily reduces via the sten- 
ographic projection to the following mixed Dirichlet-Neumann problem on a half- 
disk D := {(r,ip) E M 2 : r < 1, < tp < n} (here (r, ip) are usual planar polar 
coordinates): 

4A 

-Av = — — — — onD, v \ dlD = 0, {dv/dn)\g aD = 0. (1.3.6) 
(1 + r^y 

Here d ± D := {(r, 0) : r e (0, 1)} U {(1, ip) : - tt/2| < vr/4} and 9 2 D : = 
{(r, vr) : r G (0, 1)} U {(1, t/>) : |t/> - tt/2| > tt/4}. 

Problems (11.3.31) and (11.3.611 are quite remarkable in their own right — each 
of them is an example of a mixed Dirichlet-Neumann boundary value problem 
whose spectrum is invariant under a swap of Dirichlet and Neumann boundary 
conditions. Namely, the spectrum of (ll.3.6|) coincides with the spectrum of 

4A 

-Au = j— v or\D, v\ d2D = , (dv/dn)\ dlD = . (1.3.7) 



We refer to |JLNP| for a further discussion on Dirichlet-Neumann swap isospec- 
trality. 




Figure 1: Geometry of boundary value problems ()1.3.6|) (left) and (II .3.7)1 (right). 
Here and further on, the solid red line denotes Dirichlet boundary condition and 
the dashed blue line — the Neumann one. 



Remark 1.3.8. One can check that a surface with a finite number of conical 
singularities can be approximated by a sequence of smooth surfaces of the same 
genus and area in such a way that the corresponding sequence of the first non- 
zero eigenvalues converges to Ai on the original surface. Thus, Conjecture II. 3. II 
means that (jTTQj) is sharp in the class of smooth metrics, although the equality 
is not necessarily attained. For a general result about the convergence of the 
whole spectrum see [RoJ. 
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2 Symmetries 



2.1 Hyperelliptic involution 

Let T : V — > V, T 2 = Id be a map that intertwines the preimages of points of 
§ 2 under a two-sheeted covering II : V — > S> 2 . Clearly, the Laplace operator A 
commutes with T. 

By the spectral theorem, we can consider separately the restrictions of the 
Laplacian onto the spaces of functions which are either even or odd with respect 
to T. The even functions on V can be identified with the functions on § 2 . 
Therefore, as Ai(§ 2 ) = 2, we have Ai("P) < 2, and the equality in ()1.3.2|) will 
be achieved if and only if the first eigenvalue A° dd of the Laplacian acting on the 
odd subspace satisfies A° dd > 2. 

2.2 Isometries of V 

Consider the following isometries of 8 2 (as usual, we identify S 2 and C by sten- 
ographic projection): 

a x :z^z or (%, 77, f) (x, -77, f), 
cr 2 : z (->• -2 or (x,?7,0 >-> (~X, ?7, 0> ( 2 - 21 ) 
<r 3 :z^l/z or (x,r?,0 ^ (x^^-O- 
Here z = x+iy is a point in the equatorial plane upon which a point (%, 77, £) G S 2 
is stereographically projected. 

The hyperelliptic involution T is given by T : (2,10) — > (z,—w). For 1 < 
j < 3, a symmetry Oj of S 2 has two corresponding symmetries Sj and T o Sj 
satisfying 

n o Sj = n o t o Sj = Oj o n . (2.2.2) 

Those symmetries, with account of ()1.2.1|) are given by the explicit formulae 

si : (z, w) i-> (z, z/w), 

s 2 : (z,w) 1 — ► (— z, iw), (2.2.3) 
s 3 : (z,w) 1 ► (l/z,w/z). 

As an illustration, we demonstrate how the last of these formulae is obtained: if 
u, 2 = then by (jlTTjl . 




1\ _ 1 (l/^-e"/ 4 )(l/z-e 37ri / 4 ) _ 
1/ I(l/z + e™/ 4 )(l/z + e 3 ™/ 4 ) ~ z 2 
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thus giving the expression for s 3 . 

It easily seen that all Sj commute with T and satisfy 



Sj 2 = Id, J = 1,2,3; 
sis 3 = S3S1, s 2 s 3 = s 3 s 2 ; (2.2.4) 
s 2 si = Tsis 2 . 

Remark 2.2.5. In the proof of Theorem 11.3.51 we will use only the symmetries 
si and s 3 . Calculations for s 2 are presented for the sake of completeness (see 
Remark EH) . 

2.3 Fixed point sets of isometries 

Let Fix(S') denote a fixed point set of a mapping S. As easily seen from (I2.2.1|) . 
the sets Fix(crj), for j = 1, 2, 3, lie in the union of the coordinate lines and a unit 
circle of C, and we introduce the following notation for future reference. The 
coordinate lines are divided into two rays each by the ramification point r := 0, 
and we denote 

ai := {z = t, t > 0} , a 2 := {z = it,t > 0} , 

(13 := {z = t, t < 0} , CI4 := {z — it, t < 0} . 

The circle is divided into four arcs by the ramification points r\ := e~ m ^ 4 , 

r 2 := e 7 ™/ 4 , r 3 := e 3?n / 4 , and r 4 := e ~ 3m / 4 t and we denote the arcs by 

a k+A ■= {z = e twi/4 , t E (2k — 3,2k — 1)} , fc = 1,2,3,4, 

so that the arc a 5 goes from r 1 to r 2 , the arc a 6 goes from r 2 to r 3 , the arc a 7 
goes from r 3 to r 4 , and finally a 8 goes from r 4 to r x . 

In this notation, the fixed point sets Fix((jj) are written as 

Fix(o"i) = a± Ua 3 , Fix(a 2 ) = a 2 Ua4 , Fix(a 3 ) = a^Ua^Ua-jUa^ . (2.3.1) 

Note that each of the rays aj (j = 1, 2, 3, 4) intersects an arc a J+4 at a single 
point which we denote zf 

zi = l, z 2 = i, z 3 = -l, z 4 = -i, 

see Figure 121 
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Figure 2: Ramification points, rays, arcs and intersections 



2.4 Fixed point sets of s\, S2, S3 

Each of the points Zj has exactly two pre-i mages pf 1 ^ := (wj m \zj) G n -1 ^, 

m— 1,2, where toj 1 ' are the solutions of the equation (wj) 2 = F(zj), with F 
given in (ll.2.1|) . These solutions are easily found from ()1.2.1|) : we are of course 
at liberty to choose which of the two solutions is denoted wf and which is 
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denoted w[ 2 \ For definiteness we set 



w 



w 



(i) 



(i) _ 
1 ~~ 



y/2' 



w 



(1) 



w 



(1) 



2-V2 



(2) 



W. 



(2) 



(2) 



W 



(2) 



l + i 



2 + V2' 

-i; 

1 - i 

2- y/2' 



(2.4.1) 



For future use, we need to know the images of points pj m ^ under the sym- 
metries si, I = 1,2,3. These are easily calculated from ()2.2.3j) : it turns out that 



sip 



- with some indices k G {1, 2, 3, 4}, n G {1,2}. The results of the 



calculations are summarized in the following Tabled] 



(j, m) 


(k,n) 




I = 1 


/ = 2 


/ = 3 


(1,1) 


(1,1) 


(3,1) 


(1,2) 


(1,2) 


(1,2) 


(3,2) 


(1,1) 


(2,1) 


(4,1) 


(2,1) 


(2,1) 


(2,2) 


(4,2) 


(2,2) 


(2,2) 


(3,1) 


(3,2) 


(1,2) 


(3,2) 


(3,2) 


(3,1) 


(1,1) 


(3,1) 


(4,1) 


(2,1) 


(4,2) 


(4,1) 


(4,2) 


(2,2) 


(4,1) 


(4,2) 



Table 1: The points and their images pj^ = sip^ under symmetries 
S[. The table lists the pairs of indices (j,m) and the resulting pairs (k,n) for 
I = 1,2,3. Note also that T acts by interchanging the second indices 1^2. 



For each of the rays or arcs a^, k — 1, . . . , 8 defined in the previous section, 
its pre-image n~ 1 a/ has two connected components which we denote b^\b^ 

related by b { ™ ] = Tb { ™'\ m,m f G {1,2}, m =/ ml. Again, the choice of which 
component we denote by an upper index (1) is up to us and in order to fix 
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the notation we postulate that bjj 1 ' 3 w^T , k! = ((k — 1) mod 4) + 1, e.g. 

We now have at our disposal all the information we need in order to obtain 
the fixed point sets of s\, s 2 , S3. We start with the following two simple Lemmas. 

Lemma 2.4.2. IIFix(s i ) C Fix(a j ). 

Proof. Let z G IIFix(sj). Then there exists p eV such that lip = z and Sjp = 
p. Thus Hsjp = z and by ()2.2.2|) crjUp = UjZ = z, so that z G Fix(<x,). □ 

Lemma 2.4.3. Let C Fix(cTj). Then, form = 1,2, either b^ C Fix(sj) or 
CFix(To Si ). 

Proof. We have = Ofe, so that crjll&j^ = 0^ = a k , and so by (12.2.21) 

Ilsjl)|, m ' = afe = 116^ = nT6^ m . The result follows from the obvious observa- 
tion: if Ila = lip, then either a = (3 or a = T(5. □ 

The lemmas lead to the following 

Proposition 2.4.4. 



Fix(si) 




U6f\ 


Fix(Tsi) 


= 4 1} 


u4 2) , 


Fix(s 2 ) 


= b? 


u 4 2) , 


Fix(Ts 2 ) 


= &?> 


u bf , 


Fix(s 3 ) 


= tf } 


u 4 2) u 6« u 4 2) 


Fix(Ts 3 ) 


= tf } 


u u 4 1} u 4 2) 



Proof. By Lemmas 12.4.21 and 12.4.31 for any given j the fixed sets Fix(sj) and 
Fix(Tsj) consist only of the pre-images of the components a k of the correspond- 
ing fixed sets Fix(<7j) (given by ()2.3.1|) ). However we still need to describe which 

component b k m \ m = 1,2, lies in Fix(sj-) and which in Fix(Tsj). As each com- 
ponent b k (m) is uniquely determined by the point w k m ^ given by (12.4. ip . it is 
sufficient just to check in Table [TJ whether SjW^ = or TsjW^ = w k m \ 

For example, to find Fix(s 2 ) we need only to inspect b^ and b^™\ As, by 
Table[TJ s 2 w^ = t<4 anc ' Ts 2 w < f^ = w^ 1 ' , we have Fix(s 2 ) = 62 Ufe 2 and 
Fix(Ts 2 ) = &4 U 64 . The rest of Proposition 12.4.41 is obtained in the same 
manner. □ 
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3 Proof of Theorem 11.3.5 

We divide the proof of Theorem 11.3.51 into several steps. 



3.1 Even eigenfunctions with respect to T 

Consider the subspace V + C L 2 (V) consisting of all even eigenfunctions with 
respect to T. Any such eigenfunction has a well-defined projection on § 2 . There- 
fore, if there exists a first eigenfunction of V that belongs to V + , its projection 
is an eigenfunction on S 2 and hence the corresponding eigenvalue is greater or 
equal than two (recall that Ai(S 2 ) = 2). Hence, in this case the Conjecture [OH] 
is verified. 

3.2 Use of symmetries «i,S3. 

Denote by G 13 the subgroup of the automorphism group of V generated by the 
symmetries 

{T, si,s 3 }. 

It follows from ()2.2.4|) that G\z is commutative. Note also that all the elements 
of G 13 commute with the Laplacian on V. Therefore, we can choose a basis 
of L 2 {V) consisting of joint eigenfunctions of all s G Gi 3 and A. Given a joint 
eigenfunction / of all s G G13, we denote by s) the corresponding eigenvalue 
of s, i.e. 

f(sx) = /i(f, s)f(x) . 
Since s 2 = T 2 = Id for j = 1, 3, we see that s) = ±1 for all s G G13. 

3.3 Odd eigenfunctions with respect to T 

Consider now the space C L 2 {V) consisting of all eigenfunctions of the 
Laplacian which are odd with respect to T. Let 0i be a joint eigenfunction of 
{T, si, S3, A}, corresponding to the smallest eigenvalue of A |y_ 

Now, since fi((pi,T) = —1 and s\T = T, we have fJ>(4>u s%T) = 

/i(0i,T) = -1, and similarly s 3 )M0i, s s T ) = ~ l - 

Without loss of generality we may assume that /x(</>i,si) = — 1. We recall 
from section 12.31 that the fixed point set Fixsx consists of the arcs 6^ ,6^ . 
Thus 4>i must vanish on these arcs. 
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Consider now the symmetries s 3 , s 3 T. We must have one of the following 
two cases: 

i) s 3 T) = -1, s 3 ) = 1; 

ii) /u(0i,s 3 ) = -1, //(0i,s 3 T) = 1. 
Consider first Case i). 

Proposition 3.3.1. In Case i) the function <pi vanishes on the arcs 

6 (1) ft (2) ft (l) 6 (2) ,(1) ,(2) 

and /ts normal derivative d n (pi vanishes on the arcs 

,(1) ,(2) ,(1) ,(2) ,(1) ,(2) 
°3 > °3 > °6 ) °6 ) °8 ) °8 ' 

Proof. By Proposition 12.4.41 the fixed-point set of s 3 T consists of the arcs 
65,65,67,67. Accordingly, 0i vanishes on all those arcs, as well as on b\,b\. 
Moreover, <f>\ has /i(0i, s 3 ) = siT) = 1. It follows that the normal deriva- 
tive of d n (pi vanishes on the fixed-point sets of those symmetries. It remains to 
apply once more Proposition 12.4.41 in order to complete the proof. □ 

Consider next Case ii). 

Proposition 3.3.2. In Case ii) the function <pi vanishes on the arcs 

,(1) ,(2) ,(1) ,(2) ,(1) ,(2) 
u 3 ' °3 ' u 6 1 °6 1 U S 1 °8 ) 

and its normal derivative d n (pi vanishes on the arcs 

6 (i) 6 (2) 6 (i) 6 (2) fc (i) fc (2) 

Proposition 13.3.21 is proved in the same way as Proposition 13.3.11 

3.4 Final step of the proof 

Since (pi is an odd function with respect to the hyperel Ii ptic involution T, its 
projection upon § 2 is not well-defined. However, the projection of \(pi\ to § 2 is 
well-defined. Denote it by 

In Case i), the function tpi can be chosen as a test function for the mixed 
Dirichlet-Neumann boundary value problem (ll.3.3j) . Assume now Conjecture 
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11.3.41 is true and the first eigenvalue of (ll.3.3|) satisfies Ai > 2. Then the 
Rayleigh quotient of ipi and hence of 0i satisfies the same inequality. But this 
means that ipx cannot be the first eigenfunction on V since we get a contradiction 
with ()1.1.3|) . Therefore, the first eigenfunction of V is even with respect to T, 
and as was shown in section I3~T1 this implies Conjecture 11.3.11 

Similarly, in Case ii), the function ^ can be chosen as a test function for 
the mixed Dirichlet-Neumann boundary volume problem which is obtained from 
(11.3.3ft by swapping the Dirichlet and the Neumann conditions. However, it 
was shown in jJLNPj that this problem is isospectral to (jl.3.3p . Therefore, 
repeating the same arguments as above we prove that Conjecture 11.3.11 holds. 
This completes the proof of Theorem 11.3.51 □ 

Remark 3.4.1. In the proof of Theorem 11.3.51 we have used only the symmetries 
si and s 3 . Alternatively, we could have used s 2 and S3. One can check directly 
using Proposition I2.4.4l that applying s 2 one obtains a mixed Dirichlet-Neumann 
boundary value problem which is equivalent to (11.3.3j> and hence no additional 
information about the first eigenfunction is obtained. 

3.5 A family of extremal surfaces of genus two 

The purpose of this section is to prove the following 

Corollary 3.5.1. Coniecture \1.3.4\ implies that there exists a continuous family 
Vt of surfaces of genus 2 such that A 1 Area('P t ) = I671". 

Proof. Consider the Riemann surface V t defined by the equation 



where t G (0, 7r/2). Note that V n /4 = V. It is easy to see that for any t, Vt is 
symmetric with respect to s\ and S3. Arguing in the same way as in the proof of 
Theorem 11.3.51 and using a stereographic projection, we reduce the problem on 
Vt to the following two mixed Dirichlet-Neumann boundary value problems on 
the half-disk D: 



~ Av = T\ 4- 2^ on D ' = ' ( dv / dn )\d2(t) = • (3.5.2) 



~ Av = h j- 2\2 V on D ' v Wt) = , (dv/dn)\ dl(t) = . (3.5.3) 




4A 



and 



4A 
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Here d^t) : = {(r,0) : r E (0, 1)} U {(1, vjj) : \ip-ix/2\ < t} and d(t)D : = 
{(r, tt) : r G (0, 1)} U {(1, : tt/2 > |^ - tt/2| > *}■ 



We remark that for t ^ 7r/4 these two problems are not isospectral. Using 
Dirichlet-Neumann bracketing it is easy to see that (I3.5.2|) has a smaller first 
eigenvalue than (J3.5.3P if t < n/A and a larger one if t > 7r/4. Denote the 
minimal first eigenvalue of the two problems by Ai(i). According to Conjec- 
ture 2 and numerical calculations, Ai(7r/4) > 2. Since the first eigenvalues of 
both problems depend continuously and monotonically on parameter t, and since 
A^O) = Ai(7r/2) = 0.75 (see section [T3j), there exist numbers t\ E (0,7r/4) 
and t* 2 E (7r/4,7r/2) such that Ai(tJ) = Ai(^) = 2 and so Ai(i) > 2 for 
t E [£*,£•>]■ Arguing is above, we deduce that for all surfaces V t corresponding 
to these values of t, estimate (ll.l.3|) is sharp. This completes the proof of the 
theorem. □ 

Corollary 13.5.11 implies that 16n is a degenerate maximum for AiArea(M) 
for surfaces of genus two. This is not the case for surfaces of lower genus on 
which the metric maximizing the first eigenvalue is unique. Note also that the 
extremal metrics in genera zero and one are analytic, while the surfaces Vt have 
singular points. 

4 Numerical investigations 

4.1 Basics of the Finite Element Method 

In this section we describe the numerical experiments used to estimate the first 
eigenvalue of (11.3.70 . 

We define the space H as the closure of {v E C°°(D)\supp v n d 2 D = 0} 
with respect to the H l (D) norm. 

The variational setting for the eigenvalue problem is to find the smallest 
eigenvalue Ael, and the associated eigenvector v E H such that for all u E H, 



We use finite elements to approximate the eigenvalues and eigenfunctions of 
(14.1.10 . The general procedure we follow is: 

1. Discretize the region D using a triangular mesh T h = IJ^Tj, with a size 
of an individual triangle r E T h parameterized by h > 0. 




(4.1.1) 
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2. Introduce a finite-dimensional subspace Vh of H, consisting of finite ele- 
ment basis functions {0j}^L\ on T h ; 

3. Denote (v h , X h ) G V h x E, with ^ = 1> 2 , • • • , ^atJ*. the solution of 
the finite-dimensional generalized eigenvalue problem 



Clearly, problem ()4.1.2j) is the discrete analog of (|4.1.1|) . The eigenpair 
(vh, \h) is computed using some iterative algorithm, until a prescribed 
tolerance tol is reached. 

4. Steps 1-3 are repeated with smaller and smaller h until some other stopping 
criterion is attained. 

We now present the results of some numerical experiments based on this 
strategy. 

4.2 Conforming finite elements 

In the first set of experiments, the choice of approximating subspaces Vh was a 
sequence of Pl-conforming finite element spaces (Courant triangles), see [BrJ. 
This means that for a given h > 0, and a triangulation T h of the domain, 

Vh ■= {v G H; v\ T = polynomial of degree < 1 for every r G T h } . 

The discrete generalized eigenvalue problem for each h was solved using an 
Arnoldi iteration with shift 2.6. For details on the Arnoldi iteration, see, e.g., 



Experiment 1: The initial triangulation is based on a graded mesh, with 
more triangles located near the singularities of the eigenfunction (see Figure |3J). 

The refinement strategy was based on simply subdividing each triangle in Th 
into 4 while preserving the quality of the mesh, yielding a new triangulation %i/2- 
The eigenvalue solver was run until a tolerance of tol = 10~ 16 was achieved. The 
computation was performed using FreeFem++ for generating the finite elements 



A h v h = KB h v h 



(4.1.2) 



where 




(4.1.3) 



[^ovLl ITrBa] . 
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Figure 3: A finite element mesh 





Err 


N h 


No. of 
Triangles 


No. of 
Arnoldi iterates 


2.45590105457 


2.00434573363e-05 


169 


288 


11 


2.36301118569 


1.32592470511e-06 


625 


1152 


11 


2.32089716556 


8.16135748742e-08 


2401 


4608 


12 


2.30111238184 


4.80693483786e-09 


9409 


18432 


11 


2.29161462311 


2.79565739833e-10 


37249 


73728 


11 



Table 2: Using Pl-conforming finite elements 



and the meshes, and ARPACK for the eigenvalue solve. The meshes were refined 
until the measure of error, 



Err:= / | Vu h \ 2 - 4X h — dD 
Jd [l+r 2 Y 



satisfied |Err| < 5 x 10~ 10 . The results are tabulated in Table El 

Experiment 2: This experiment was conducted using MATLAB's finite el- 
ement package PDEToolbox, and the eigenvalue solve was performed using 
ARPACK routines. A sequence of triangular meshes was created, starting from 
the coarsest mesh, and refining 5 times. The major difference between this and 
the previous experiment is in the manner in which the zero Dirichlet data is 
enforced. 
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Nu 


No. of 
Triangles 


2.55310562723060 


77 


126 


2.40400118356918 


279 


504 


2.33742285062686 


1061 


2016 


2.30582934149898 


4137 


8064 


2.29039772121374 


16337 


32256 


2.28276090970583 


64929 


129024 


2.27895954902635 


258881 


516096 



Table 3: Using Pl-conforming finite elements in MATLAB 



4.3 Nonconforming finite elements 

In the second set of experiments, we used Pl-noncon forming finite elements 
(Crouzeix elements), see [Br]. These are defined as 

V h := {v G L 2 (D);v\ T is linear for each r G T h , 

v is continuous at the midpoints of triangle edges} 

for a given h > and a triangulation T h . Note that Vh is not a subspace of H; for 
more information on the use of nonconforming elements in eigenvalue problems, 
see jArDu) . As before, the discrete generalized eigenvalue problem is solved 
using an Arnoldi iteration with a shift of 2.2 until a tolerance of tol = 10~ 16 
is achieved. The refinement strategy was to subdivide each triangle in Th into 4 
subtriangles, yielding a new mesh T h /2- The meshes were refined until a measure 
of error 2 

Err:= / | Vv h \ 2 - 4A h — dD , 

Jd {i+r 2 y 

satisfied |Err| < 5 x 10~ 10 . The results are presented in Table |U 

In each of the experiments above, we found that the computed eigenvalues 
appeared to converge to a value greater than 2.27. The associated eigenfunctions 
also appear to converge to a function whose contour lines are shown in Figure El 
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Err 


N h 


No. of 


No. of 








Triangles 


Arnoldi iterates 


2.13042989031 


-1.5494060025e-05 


169 


288 


11 


2.20743747322 


-7.89999667122e-07 


625 


1152 


11 


2.24561396752 


-3.8455570927e-08 


2401 


4608 


11 


2.26440630518 


-1.87263030138e-09 


9409 


18432 


11 


2.27364314423 


-8.78059287464e-ll 


37249 


73728 


11 



Table 4: Using Pl-nonconforming finite elements 




Figure 4: Contour lines of the first eigenfunction 
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